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ABSTRACT 


Five  synoptic  space  sections  along  158  W  longitude  between  Hawaii 
and  the  Aleutian  Islands  were  developed  from  data  collected  by  air- 
borne expendable  bathythermographs  during  experiment  PARKA,  a  research 
project  sponsored  by  the  U.  S.  Navy  in  1968.   The  sections  are  examined 
for  spacial  and  temporal  variation  in  thermal  structure  and  geostrophic 
surface  velocity. 

Two  recently  developed  analysis  techniques  are  employed.   Denner's 
T-S  gradient  method,  wherein  thermal  and  haline  contributions  to  total 
geostrophic  velocity  are  distinguishable,  expedites  calculations  and 
results  in  velocity  fields  comparable  to  those  developed  by  the  dynamic 
method.   Thermocline  parameters  are  developed  using  Boston's  objective 
definition  of  the  thermocline,  a  statistical  curve-fitting  technique 
which  develops  the  notion  of  a  Gaussian  thermocline. 

Gross  features  of  thermal  structure  remain  fairly  consistent 
during  the  heating  season;  however,  thermal  fronts  are  observed  to 
vary  in  time  and  space.   The  distribution  of  isothermal  lines  with 
latitude  suggests  the  possibility  of  a  Taylor-column  effect  slightly 
north  of  Hawaii. 
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I.   INTRODUCTION 

The  genesis  of  this  investigation  was  experiment  PARKA  -  a  research 
project  sponsored  by  the  Oceanographer  of  the  Navy  to  investigate  the 
acoustic  environment  in  the  mid  North  Pacific  Ocean.   The  objectives  of 
PARKA  were  many  and  varied  and  are  Largely  inappropriate  for  discussion 
here.   More  to  the  point  is  the  large  body  of  environmental  data  accumu- 
lated during  the  execution  of  the  project. 

The  principal  phases  of  PARKA  were  conducted  in  August  and  September 
1968.   During  this  period  more  than  2,000  observations  were  recorded 
along  a  line  approximately  corresponding  to  158  W  longitude  between 
22-55  N  latitude.   More  than  ten  surface  ships  (including  the  floating 
instrument  platform  FLIP)  together  with  Navy  patrol  aircraft  partici- 
pated in  the  experiment.   Because  of  the  acoustic  nature  of  the  mission, 
the  program  of  environmental  data  collection  was  designed  primarily  to 
determine  thermal  structure.   To  this  end,  a  variety  of  instruments  were 
used:   mechanical  bathythermographs  (BT) ,  shipborne  (XBT)  and  airborne 
(AXBT)  expendable  bathythermographs,  airborne  radiation  thermometers 
(ART),  and  in  situ  salinity,  temperature,  depth  recording  devices  (STD) . 
Salinity  structure  along  the  track  was  defined  from  climatological 
atlases  augmented  by  a  limited  number  of  real-time  profiles  recorded 
by  STD. 

Considered  in  total,  the  PARKA  data  define  the  thermal  structure 
in  the  upper  300m  of  a  hydrographic  section  extending  some  2,000  miles. 
The  section  is  unique  in  two  respects.   First  is  the  sheer  number  of 
observations  wherein  thermal  parameters  were  determined  independently 


and  often  simultaneously  by  several  different  instruments.   Secondly, 
the  sections  developed  from  AXBT  records  are  as  nearly  synoptic  as 
possible  -  roughly  70  profiles  along  the  entire  track  completed  in 
eight  hours. 

The  most  intriguing  aspect  of  PARKA  for  the  author  was  the  extensive 
use  of  the  AXBT.   The  AXBT  flights  not  only  initiated  the  instrument  in 
a  major  scientific  undertaking  but  also  produced  the  only  known  space 
sections  of  thermal  structure  on  so  large  a  scale.   The  entire  track 
between  Hawaii  and  the  Aleutian  Islands  was  traversed  five  times  -  four 
flights  during  the  actual  execution  of  PARKA  and  a  follow-on  mission  two 
months  later.   The  fifth  transit  was  completed  23  November  1968;  thus 
the  time  frame  of  the  flights  includes  the  end  of  a  heating  season  and 
the  first  one-third  of  a  cooling  season.   Five  sections  developed  from 
data  collected  during  these  flights  form  the  basis  for  this  study. 

Each  section  was  examined  to  determine  latitudinal  variation  of 
selected  isotherms  and  of  the  depths  to  the  top,  center,  and  bottom 
of  the  thermocline.   One  section  was  selected  for  geostrophic  com- 
putations by  both  the  dynamic  method  (Proudman,  1953)  and  the  T-S 
gradient  method  (Denner,  1969)  to  compare  resulting  surface  velocity 
fields.   A  correlation  between  total  surface  velocity  and  its  thermal 
contribution,  as  calculated  by  the  T-S  gradient  method,  was  attempted 
in  order  to  develop  geostrophic  velocity  fields  in  the  remaining 
sections  based  solely  on  thermal  structure. 

With  the  spacial  variation  of  the  indicated  parameters  established, 
the  five  sections  were  considered  together  in  order  to  develop  the 
objectives  of  this  thesis.   Briefly  stated,  these  objectives  are:   to 
examine  the  spacial  and  temporal  variation  of  the  top,  center  and 
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bottom  of  the  thermocline ;  to  examine  the  spacial  and  temporal  variation 
of  selected  isotherms,  correlating  the  results  to  thermocline  variations; 
to  test  the  T-S  gradient  method  as  a  device  for  rapid  and  realistic 
dynamic  calculations;  and  to  examine  the  AXBT  as  a  scientific  instrument. 

In  developing  these  objectives,  two  fairly  recent  inovations  are 
employed.   The  first,  already  mentioned,  is  the  T-S  gradient  method 
wherein  the  thermal  and  haline  contributions  to  geostrophic  velocity 
are  distinguishable.   The  second  is  the  notion  of  a  "Gaussian  thermo- 
cline" developed  by  Boston  (1966).   In  order  to  bring  the  results 
developed  from  these  two  techniques  into  sharper  focus  -  as  well  as 
to  evaluate  the  reliability  of  the  AXBT  -  this  thesis  considers 
exclusively  the  temperature  and  depth  data  provided  by  airborne  thermo- 
grams.  The  assumption  is  that  if  reasonable  results  can  be  developed 
solely  from  temperature  data,  the  subjective  application  of  atmospheric 
effects  at  the  air/sea  interface  will  only  enhance  the  value  of  the 
results  . 
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II.   DISCUSSION  OF  THE  DATA 

Temperature  data  were  collected  during  five  flights  by  Navy  P-3 
patrol  aircraft.   Altitude  and  airspeed  varied  somewhat  among  missions, 
depending  on  flying  conditions,  but  remained  well  within  limits  for 
efficient  operation  of  the  AXBT.   Altitudes  generally  ranged  between 
2,000-5,000  feet  and  airspeeds  between  200-250  knots. 

A.   AREA  COVERAGE 

Each  of  the  five  flights  originated  in  the  Aleutian  Islands  and 
terminated  in  Hawaii.   Flight  plans  were  the  same  for  each  mission. 
The  track  extended  1,985  nautical  miles  along  157  50'W  longitude 
between  55  05'N  and  22  00'N  latitude.   Figure  1  is  a  representation 
of  the  cross  section  along  the  track.   Take-off  time  was  scheduled  to 
ensure  arrival  of  the  aircraft  on  the  track  one  hour  either  side  of 
0800  local  time.   Optimum  spacing  between  observations  was  set  at 
25nmi  which  corresponds  to  an  instrument  launching  about  every  six 
minutes . 

All  flights  conformed  to  the  track  within  navigational  error  and 
the  instrument  launching  schedule  was  followed  as  closely  as  practi- 
cable.  Spacing  between  temperature  profiles  along  the  track  was 
ultimately  a  function  of  the  quality  of  the  thermograms;  only  those 
traces  extending  to  300m  or  greater  were  considered.   The  lack  of 
uniformly  spaced  observations  does  not  inhibit  development  of  the 
objectives  outlined  earlier. 

Table  I  is  a  summary  of  pertinent  time  and  space  factors  describing 
each  section;  the  indicated  number  of  observations  are  those  that  met 
the  300m  criterion  and  were  free  from  obvious  errors. 
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SECT 

DATE 

TIMES 
(local) 

COVERAGE 
(°N  lat) 

OBS 

1 

19  Aug 

68 

0809-1630 

54.1-22.3 

70 

2 

22  Aug 

68 

0807-1621 

54.7-22.6 

72 

3 

2  Sept 

68 

0810-1636 

54.6-22.1 

67 

4 

4  Sept 

68 

0810-1649 

54.7-22.2 

69 

5 

23  Nov 

68 

0840-1746 

53.9-22.0 

55 

TABLE  I.   Time  and  space  dimensions  of  synoptic  sections 
along  158°W  longitude. 


B.   INSTRUMENTATION 

The  instrument  package  carried  aboard  the  aircraft  included 
AN/SSQ-36  Bathythermograph  Transmitter  Sets  (herein  referred  to  as 
AXBT)  and  an  Anadex  Converter-Rustrack  Recorder  Assembly  (NADC  Con- 
verter) used  in  conjunction  with  an  AN/ARR-52  Radio  Receiving  Set. 

As  the  name  implies,  the  AXBT  is  an  expendable  aircraft  launched 
device  designed  to  provide  a  sea  temperature  profile  from  the  surface 
to  a  depth  of  1,050  ft.   Major  components  of  the  instrument  are  a 
temperature  probe,  an  audio  oscillator  controlled  by  the  probe,  a 
timing  circuit,  batteries  and  a  VHF  transmitter.   The  total  assembly 
is  housed  in  a  standard  size  sonobuoy  case  and  weighs  18  pounds. 

The  AXBT  may  be  launched  from  virtually  any  altitude;  however, 
something  less  than  10,000  ft.  produces  best  results.   At  an  altitude 
of  3,000  ft.,  for  example,  a  complete  thermogram  is  recorded  aboard  the 
aircraft  about  5.5  minutes  after  the  instrument  is  launched;  at  1,000  ft 
the  interval  would  be  about  4.7  minutes.   The  instrument  is  stabilized 
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by  a  rotochute  during  decent.   On  impact,  the  bottom  plate  is  released 
and  the  sea -water  batteries  become  active.   Approximately  75  seconds 
later,  the  timing  circuit  releases  the  temperature  probe  which  sinks  at 
a  uniform  rate  of  five  feet  per  second.   An  audio  signal  whose  frequency 
is  proportional  to  temperature  modulates  the  VHF  carrier  frequency  during 
descent  of  the  probe.   The  audio  frequency  varies  from  1,300  to  2,700Hz 
corresponding  to  a  temperature  range  of  25  to  95  F.   Primary  voltage  is 
removed  from  the  transmitter  about  six  minutes  after  entry  into  the  water; 
a  water-soluable  plug  dissolves  soon  thereafter  allowing  the  instrument 
to  sink. 

The  Anadex  Converter  receives  the  demodulated  signal  directly  from 
the  ARR-52  radio  receiver  and  converts  it  to  a  voltage  proportional  to 
the  frequency  of  the  incoming  signal.   The  converter  output  signal  is 
supplied  to  the  Rustrack  Recorder  where  the  signal  amplitude  developed 
by  the  converter  is  visually  displayed  on  a  recording  chart.   The  chart 
drive  is  set  at  three  inches  per  minute  to  allow  in-flight  interpretation 
if  required. 

The  effectiveness  of  the  AXBT  is  discussed  in  section  V. 

C.   PROCESSING  AND  ANALYSIS 

At  the  completion  of  each  flight,  the  thermograms  were  delivered  to 
a  representative  of  Fleet  Weather  Center,  Hawaii.   All  traces  are  for- 
warded eventually  to  the  Fleet  Numerical  Weather  Central  (FNWC)  in 
Monterey,  California,  where  they  are  digitized  and  stored  on  tape. 
The  ultimate  repository  for  the  traces  is  the  Naval  Oceanographic 
Office,  Washington,  D.  C. 
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This  processing  has  been  completed  for  the  thermograms  defining 
sections  1-4  (see  Tab.  I).   Traces  from  section  five  were  interpreted 
at  the  Hawaii  Institute  of  Geophysics  and  results  transmitted  electron- 
ically to  FNWC.   Digitized  versions  of  these  traces  were  not  available 
at  this  writing. 

All  observations  are  coded  as  a  series  of  linear  segments.   In  the 
case  of  the  digitized  data,  the  CDC  6500  computor  at  FNWC  was  programed 
to  report  temperature/depth  pairs  so  that  temperature  deviated  from 
linearity  by  no  more  than  0.075  C  between  reported  levels.   The  data 
transmitted  electronically  from  Hawaii  reported  depths  and  associated 
temperatures  of  observed  points  of  inflection  in  the  traces. 

All  data  were  punched  on  cards  for  processing  by  the  Naval  Post- 
graduate School  IBM  360  computor.  A  linear  interpolation  scheme  was 
devised  to  permit  temperature  determinations  at  more  convenient  depths. 

Thermocline  calculations  were  made  using  subroutines  developed  by 
Grosfils  (1968)  following  Boston  (1966).   Algorithms  for  computing 
geostrophic  velocity  were  written  by  Mr.  R.  L.  Middleburg  of  the 
Department  of  Oceanography,  Naval  Postgraduate  School,  and  are  pre- 
sented by  Denner  (1969).   The  contours  of  isothermal  lines  were  plotted 
by  Subroutine  METMAP  from  the  IBM  library  at  the  Naval  Postgraduate 
School. 
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III.   GEOSTROPHIC  MOTION  ALONG  158°W  LONGITUDE 


Geostrophic  calculations  by  the  classical  dynamic  method  are 
described  in  any  of  several  textbooks  on  physical  oceanography,  for 
example  Proudman  (1953)  or  Neumann  and  Pierson  (1966).   The  limitations 
and  computational  errors  inherent  in  the  dynamic  method  are  summarized 
by  Denner  (1969).   As  an  alternative  to  the  dynamic  method,  Denner  has 
developed  a  simplified  scheme  in  which  geostrophic  velocity  is  expressed 
as  three  terms  involving  the  observed  variables,  salinity  and  tempera- 
ture.  A  useful  feature  of  the  method  is  that  the  thermal  (V  )  and 
haline  (V,  )  contributions  to  total  geostrophic  velocity  (V  )  are 
calculated  separately.   The  import  of  this  valuable  feature  is  the 
possibility  for  rapid  geostrophic  calculations  based  solely  on  thermal 

structure  -  provided  a  correlation  between  V  and  V  can  be  established. 

g      t 

Such  a  correlation  was  determined  by  Moynihan  (1968)  in  two  water  masses 
in  the  Grand  Banks  area.   He  pointed  out  the  tangible  benefits  that 
would  accrue  to  the  International  Ice  Patrol  if  geostrophic  velocity 
fields  could  be  determined  from  a  network  of  stations  occupied  by 
bathythermographs  alone,  making  use  of  the  high  ship  speeds  permitted 
by  use  of  the  XBT.   Granting  this,  even  greater  savings  in  time  and 
effort  would  result  from  the  use  of  airborne  rather  than  surface 
instrument  platforms. 

Accordingly,  this  section  briefly  describes  the  development  of  the 
T-S  gradient  method  and  applies  the  technique  to  the  AXBT  data. 

A.   THE  T-S  GRADIENT  METHOD 

Assuming  that  sea  water  can  be  considered  a  binary  fluid  system, 
Denner  followed  Reid  (1959)  in  observing  that  the  He lland-Hansen 
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equation  for  geostrophic  flow  between  two  stations  can  be  expressed: 

dP  (1) 


w~-i  [[<&,[,%+(&„ 


dS 
dnJ 


Pl 


where 


n  =  the  horizontal  distance  between  stations; 

f  =  the  Coriolis  parameter; 

a  =  specific  volume  as  a  function  of  salinity,  temperature, 
and  pressure  ; 

(V-j-V^)  =  magnitude  of  the  relative  current  normal  to  n. 
The  subscript  P  indicates  the  bracketed  quantity  is  evaluated  at  con- 
stant pressure. 

Equation  1  expresses  the  geostrophic  velocity  in  terms  of  three 
essential  elements: 

Horizontal  gradients  of  temperature  and  salinity; 

A  quantity  specifying  the  dependence  of  specific  volume  on 
temperature  at  a  given  salinity  and  pressure;  and 

A  quantity  specifying  the  dependence  of  specific  volume  on 
salinity  at  a  given  temperature  and  pressure. 

Note  that  the  specific  volume  terms  in  (1)  become  the  coefficients  of 
thermal  compressibility  and  saline  contraction  if  each  is  divided  by 
specific  volume;  compressibility  of  sea  water  enters  the  equation  in- 
directly through  the  dependence  of  these  quantities  on  pressure. 

Using  (1)  as  a  departure  point,  Denner  examined  several  possibilities 
for  developing  tractable  expressions  for  the  specific  volume  terms.   The 
approach  selected  was  to  compute  specific  volume  at  varying  values  of 
temperature,  salinity  and  pressure  using  the  Ekman  equation  of  state 
for  sea  water  and  fitting  the  results  with  polynomials  in  the  least 
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squares  sense.   The  Ekman  equation  was  selected  because  of  its  wide 

usage  rather  than  any  intrinsic  superiority  over  other  equations  of 

state . 

Values    of   specific   volume  were    computed    for   all   combinations    of 

variables   at    prescribed    intervals    between   0-30  C    temperature,    30-40   /oo 

salinity,    and   0-200bar    pressure.      The   resulting   polynomials   were   of 

the    form: 

a(T)  =  a    (S,P)   +  A(S,P)T  +  B(S,P)T2  (2) 

o  ,F  o 

a(S)T,P  =  ao(T>P)  "  C<T>P)S  <3> 

where  a   is  the  specific  volume  at  atmospheric  pressure. 

Differentiation  of  (2)  and  (3)  with  respect  to  temperature  and 
salinity  respectively  yields  expressions  appropriate  for  substitution 
into  (1).   Effecting  this  substitution  produces  the  following  simpli- 
fied expression  for  computing  geostrophic  currents,  defined  by  Denner 
as  the  T-S  gradient  method: 

P2 

(W  =  f  /  [A<s'p>  £  +  ^os.m  £  -  c(t,p)  ||]  dp         (4) 
pi  P 

He  pointed  out  that  the  interaction  between  the  thermal  and  haline  con- 
tributions to  total  flow  takes  place  through  the  coefficients  A(S,P), 
B(S,P),  and  C(T,P);  therefore,  the  first  two  terms  of  (4)  can  be  con- 
sidered the  thermal  contribution  and  the  last  term  the  haline  contri- 
bution to  total  flow.  Interpolated  values  of  the  three  coefficients  as 
functions  of  their  independent  variables  are  displayed  in  graphic  and 
tabular  form  in  Denner' s  paper. 
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1.   Computational  Form 

To  extend  application  of  the  T-S  gradient  method,  Denner 
demonstrated  that  weighted  mean  values  of  the  three  coefficients  in 
(4)  could  be  determined  for  large  oceanic  areas  to  facilitate  calcu- 
lation of  geostrophic  surface  currents  relative  to  a  selected  isobaric 
level.   Equation  4  could  then  be  rewritten  in  terms  of  finite  differences 
for  computation: 


<VVi  "  is 

where : 


K.    S   AT.AP.   +  K»   S  T.AT.AP      -   k     S   AS    AP 
-    1   i        ii  2   i     l     l      i  3    i        iiJ 


(5) 


K1 ,  K2,  K„  =  weighted  mean  values  of  A(S,P),  2B(S,P),  and 

C(T,P)  respectively  between  the  surface  and  the 
selected  reference  level  (cgs  units); 

AT.  =  average  horizontal  temperature  difference  in  the 
i-th  pressure  interval  (°C) ; 

T.  =  average  temperature  in  the  i-th  pressure  interval 
1   <°C); 

AS.  =  average  horizontal  salinity  difference  in  the 
i-th  pressure  interval  (o/oo) ; 

2 
AP.  =  the  i-th  pressure  interval  (dynes/cm  ) ; 

An  =  horizontal  separation  of  the  stations  (cm). 

Denner  assigned  fixed  values  to  K..  ,  K  ,   and  K„  relative  to  the 

l,000db  surface  in  the  Pacific  Ocean  after  observing  that  average  values 

of  A(SSP)5  B(SjP)  and  C(T,P)  calculated  for  summer  and  winter  conditions 

did  not  vary  greatly  over  the  entire  area  or  between  seasons.   The 

assigned  values  are: 

—  f\ 

Kj  =  81.0  x  10"   (cgs  units) 

K_  =  86.5  x  10"   (cgs  units) 

1L  =  -74.0  x  10"   (cgs  units) 
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Calculations  using  these  fixed  values  introduced  errors  of  not  more 
than  five  percent  in  surface  velocities  computed  from  actual  observations 
2.   Application  to  AXBT  Data 

A  preliminary  step  to  the  task  implied  in  the  heading  of  this 
section  was  to  validate  the  correspondence  between  surface  velocities 
determined  by  the  dynamic  method  and  by  the  T-S  gradient  method.   AXBT 
section  2  was  selected  arbitrarily  for  the  comparison.   Salinity  pro- 
files corresponding  to  those  of  temperature  were  supplied  by  FNWC  from 
historical  data  in  their  tape  files.   Absolute  values  of  salinity  were 
not  important  in  that  the  same  values  were  used  in  both  methods;  more 
germane  was  the  characteristic  variety  of  structural  features  provided 
by  the  profiles.   Surface  velocities  were  calculated  by  both  methods 
for  72  station  pairs ;  the  pairs  were  formed  by  coupling  adjacent  stations 
in  succession  along  the  track.   The  relative  error  between  results  for 
each  station  pair  was  calculated  by  dividing  the  difference  between  the 
two  velocities  by  the  magnitude  of  the  "dynamic"  velocity  and  expressing 
the  results  as  a  percentage.   Figure  2  illustrates  the  correlation  bet- 
ween results. 

The  agreement  was  excellent.  Over  half  of  the  73  comparisons 
differed  by  less  than  two  percent.  Relative  error  exceeded  five  percent 
only  at  those  points  where  calculated  velocities  were  on  the  order  of 
one  centimeter  per  second.  The  introduction  of  large  relative  errors 
in  regions  of  weak  currents  does  not  detract  from  the  practicality  of 
the  T-S  gradient  method.  Quite  often  it  is  sufficient  merely  to  know 
that  geostrophic  currents  are  less  than  a  few  centimeters  per  second; 
this  fact  can  be  determined  even  if  the  relative  error  is  large. 
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The  question  to  be  addressed  next  is  whether  or  not  a  definitive 

relation  exists  between  the  geostrophic  surface  velocity  and  its  thermal 

contribution.   Throughout  the  remainder  of  this  paper,  references  to 

geostrophic  velocities  or  currents  imply  values  calculated  by  the  T-S 

gradient  method. 

a.   Correlation  Between  V  and  V 

g      t 

Again,  section  2  was  selected  to  demonstrate  correlation. 

V  was  calculated  using  (5) ;  V  was  calculated  using  (5)  less  the  third 

term  on  the  right  side.   The  reference  level  was  set  at  300m  and  the 

fixed  values  of  the  K-constants  determined  by  Denner  were  used. 

Initially,  station  pairs  were  formed  in  the  same  manner  as  described 

earlier,  i.e.  successive  coupling  of  adjacent  stations.   The  calculated 

velocities  were  fitted  with  a  first-degree  polynomial  in  the  least 

square  sense  by  digital  computor  using  the  IBM  library  subroutine 

LSQPOL,,   The  resulting  linear  expression, 

V  =  -0.880  +  1.040V 
g  t 

has  a  standard  deviation  of  3.44cm/sec.   The  average  spacing  between 
stations  was  on  the  order  of  50km. 

Reed  and  Laird  (1966)  point  out  that  there  is  an  optimum 
spacing  between  stations  that  should  be  considered  when  making  geo- 
strophic calculations.   If  calculations  are  made  by  the  dynamic  method, 
the  spacing  is  determined  by  the  maximum  error  acceptable  in  the  results 
and  the  measurement  error  in  determining  the  geopotential  anomaly 
difference  of  the  two  stations.   If  the  allowable  error  in  the  results 
is  set  at  20$,  for  example,  the  anomaly  difference  should  be  at  least 
five  times  the  measurement  error.   Thus,  if  measurement  error  is  +  1.1 
dynamic  centimeters,  as  they  assumed,  stations  should  be  chosen  such 
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that  the  anomaly  difference  is  at  least  5.5  dynamic  centimeters.   Table 
II  presents  the  minimum  station  spacing  recommended  by  Reed  and  Laird 
along  two  sections  of  the  North  Pacific,  one  of  which  (160  W)  closely 
corresponds  to  the  track  considered  here. 


CURRENT 

OPTIMUM  STATION 

SPACING  (km) 

along  approx. 

along  approx. 

160°W 

175°W 

Alaska  Stream 

16 

8 

Subarctic 

183 

110 

West  Wind  Drift 

61 

39 

Gyre 

>  137 

>  137 

Kuroshio  Extension 

55 

34 

TABLE  II.   Computed  minimum  distance  between  stations 

for  near-surface  geostrophic  flow  determinations. 
(After  Reed  and  Laird,  1966.) 

The  relevance  of  the  values  in  Table  II  to  calculations 

other  than  by  the  dynamic  method  is  not  clear;  unfortunately,  time 

constraints  prevented  a  more  rigorous  investigation  of  this  point.   In 

view  of  the  objective  of  this  section  -  to  develop  a  rapid,  simple 

method  for  determining  geostrophic  velocity  fields  -  a  compromise  was 

struck.   The  computor  was  programed  to  select  stations  such  that 

spacing  was  always  greater  than  100km.   For  section  2  the  result  was 

25  station  pairs  with  an  average  spacing  on  the  order  of  140km.   The 

calculated  velocities  were  again  fitted  to  a  first-degree  polynomial. 

The  resulting  expression,  equation  (6),  has  a  standard  deviation  of 

2 .473cm/sec . 
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V  =  -0.788  +  0.997V  (6) 

g  t 

The  improved  value  of  standard  deviation  supports  the  notion  of  optimum 
spacing;  additionally,  the  goodness  of  fit  would  seem  to  suggest  that 
the  horizontal  distribution  of  temperature  provides  the  major  contri- 
bution to  geostrophic  velocity  along  158  W  longitude. 

Geostrophic  velocities  were  calculated  for  the  five  sections 
using  (6);  the  results  are  presented  in  Figs.  3-7.   Reed  and  Laird 
(1966)  made  similar  calculations  (using  the  dynamic  method)  from  data 
collected  between  1961-64  along  approximately  the  same  track.   Figure 
8  presents  their  results  compared  to  actual  current  measurements. 

B.   LATITUDINAL  VARIATION  OF  GEOSTROPHIC  FLOW 

A  general  pattern  of  eastward  flow  is  evident  in  Figs.  3-7;  rela- 
tively strong  flow  recurs  at  approximately  24,  27,  and  34-36  N  latitude. 
In  Fig.  8,  Reed  and  Laird  indicate  comparable  movement  at  24,  30  and 
39  N  latitude.   Considering  the  difference  in  data  and  analysis  technique, 
the  agreement  is  remarkably  good,  both  in  location  and  magnitude.   Two 
additional  features  which  recur  consistently  in  Figs.  3-7  are  indications 
of  weak  westerly  flow  in  the  far  north  and  stronger  westerly  flow  in  the 
far  south.   The  former  is  quite  likely  a  remnant  of  the  Alaska  Gyre, 
inducing  a  westerly  velocity  component  at  the  convergence  of  the  Alaska 
Stream  and  the  Subarctic  current  (Tabata,  1965,  Fig.  1).   The  mechanism 
producing  the  latter,  however,  is  not  clear.   This  point  will  be 
addressed  further  in  the  section  discussing  isotherm  patterns.   Aside 
from  these  features,  the  five  velocity  fields  seem  to  differ  most  with 
respect  to  season. 
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Figure  3.     Latitudinal  variation  of  geostrophLc  flow -10  AUG   1968 
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Figure  5.     Latitudinal  variation  of  geostrophic  flow  -  2  SEPT  1968 
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Figure  6.     Latitudinal  variation  of  geostrophic  flow -4  SEPT  1968 
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Figure  7.      Latitudinal  variation  of  geo   trophic  flow  -  23  NOV  1968 
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Figure  8.      Comparison  of  calculated  geostrophic  surface  currents 

versus  direct  measurements  (After  Reed  and  Laird,   1966). 
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1 .  In  the  Heating  Season 

Figures  3-6  represent  the  surface  currents  across  the  track 
during  the  final  weeks  of  the  heating  season.   In  the  mid-latitudes, 
between  33-52  N,  the  flow  appears  generally  weak  and  to  the  east.   The 
strongest  currents  occur  below  33  N  and  vary  in  direction. 

The  overall  pattern  is  dominated  by  a  strong  easterly  flow 
centered  at  27.7  N.   The  location  of  this  current  remains  fixed;  the 
magnitude  averages  about  14cm/sec  in  August  and  20cm/sec  in  September. 
Immediately  to  the  south,  the  flow  pattern  alternates  to  the  west  then 
to  the  east.   The  westerly  current  is  strongest  on  19  August  (Fig.  3) 
at  15cm/sec  located  at  25.4  N.   During  the  remainder  of  the  heating 
season,  the  location  varies  slightly  but  the  speed  reduces  to  an 
average  of  about  five  cm/sec.   The  easterly  flow  at  24  N  averages 
approximately  14cm/sec  in  August  and  six  cm/sec  in  September. 

The  most  variable  feature  in  Figs.  3-6  is  the  pattern  of 
westerly  flow  just  north  of  29  N.   On  19  August  and  2  September  (Figs. 
2  and  4)  there  are  westerly  toungues  at  30  N  (average  speed  7.3cm/sec) 
which  vanish  in  the  succeeding  two  and  three  days.   Additionally,  the 
westerly  currents  located  at  32.4  N  in  August  shift  one  degree  to  the 
south  by  September. 

A  final  feature  of  interest  is  the  sudden  appearance  of  a 
relatively  strong  eastward  current  at  44  N  on  4  September  (Fig.  6). 
This  current  is  marginally  present  in  previous  figures,  and  either 
disappears  or  moves  south  by  November. 

2 .  In  the  Cooling  Season 

Figure  7  presents  the  velocity  field  after  about  two  months 
of  cooling.   Not  unexpectedly,  the  pattern  has  changed  considerably. 
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The  only  significant  westward  flow  is  in  the  extreme  south;  the  weak 
westward  flow  persists  at  50  N.   Along  the  remainder  of  the  track, 
eastward  movement  has  intensified,  developing  fairly  strong  currents 
at  24,  29,  37,  and  41  N  latitude.   This  corresponds  to  the  general 
eastward  drift  characteristic  of  the  central  subarctic  Pacific. 
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IV.   THERMAL  STRUCTURE  ALONG  158°W  LONGITUDE 


Thermal  structure  can  be  defined  as  the  relationship  between 
temperature  and  depth  as  a  function  of  time  and  space.   In  this  con- 
text, the  five  AXBT  sections  are  particularly  adaptable  to  an  investi- 
gation of  thermal  structure.   There  are,  of  course,  several  descriptive 
parameters  and  analytic  techniques  available  to  describe  structure;  two 
have  been  chosen  to  present  the  AXBT  data.   The  first  and  most  straight- 
forward is  to  plot  selected  isopleths  of  temperature  (isotherms)  along 
each  section.   The  second  is  to  show  the  spacial  variation  of  the  thermo- 
cline  by  plotting  calculated  depths  to  the  top  (Z  ),  center  (Z  )  and 
bottom  (Z,  )  of  the  thermocline  versus  latitude. 

D 

The   definition  and   determination   of   the  Z-parameters    follows    Boston 

(1966).      Boston's   method    of   analysis    is    essentially   a   curve-fitting 

technique   making   use    of   some    fundamental   statistical    theorems.      He 

introduced    the    notion   of  a    "Gaussian  Thermocline",    developed   expressions 

for   associated    parameters,    and    illustrated    the   application   of    these 

parameters    to   non-Gaussian  distributions. 

Z      and   Z,  ,    defined    in   terms    of  Z    ,    approximate    the   depths    to    the 
t  b  c 

upper   and    lower    limits    of   a    layer   of  water  with  a    large    temperature 

gradient.      The    fundamental   parameter    is   Z    ,      To  define   Z      it    is 

r  c  c 

necessary  to  introduce  a  function  Z(T)  which  describes  a  distribution 
of  depths  as  a  function  of  temperature.   This  distribution  is  developed 
from  original  data  by  forming  an  ordered  set  of  depths,  beginning  at  the 
surface,  corresponding  to  temperatures  separated  by  a  constant  increment, 
Z   is  then  defined  as  the  depth  corresponding  to  the  point  of  inflection 
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in  the  distribution  function  Z(T).   Details  of  the  method  are  summarized 
briefly  in  the  following  sections. 

A.   THE  GAUSSIAN  THERMOCLINE 

A  Gaussian  thermocline  is  defined  by  a  distribution  of  temperature 
as  a  function  of  depth,  T(Z),  whose  frequency  function  (i.e.  whose 
first  derivative)  forms  a  normal  or  Gaussian  curve.   Characteristic  of 
a  Gaussian  curve  is  symmetry  about  the  first  moment  and  symmetrically 
located  points  of  inflection  and  points  of  maximum  rate  of  change. 
Making  use  of  these  features,  Boston  demonstrated  that  Z   corresponds 
to  the  first  moment  (Z)  and  that  Z   and  Z,  correspond  to  points  of 
maximum  rate  of  change  (Z  +  v3S),  where  S  is  standard  deviation. 
Implicit  in  this  development  is  a  normalized  temperature  field  in  which 
successively  decreasing  values  of  temperature  are  associated  uniquely 
with  increasing  depths  (i.e.  the  distribution  functions  T(Z)  and  Z(T) 
are  coincident). 

Grosfils  (1968)  has  developed  numerical  methods  for  analyzing  a 
Gaussian  thermocline  by  digital  computor.   The  requisite  inputs  are 
the  surface  and  bottom  temperatures  and  depths  read  at  a  constant 
temperature  interval.   Basically,  the  technique  is  to  compute  the  first 

and  second  moments  of  the  frequency  function  which  define  mean  (Z)  and 

2 
variance  (S  )  respectively.   He  points  out,  however,  that  the  error 

inherent  in  transforming  from  integrals  to  summations  to  compute  moments 
is  closely  related  to  the  temperature  increment  chosen  to  define  the 
depth  vector.   The  error  is  most  apparent  in  the  value  of  standard 
deviation  (S)  developed  from  the  second  moment;  therefore,  he  intro- 
duced a  weight  factor  K  which  varies  as  the  temperature  increment.   This 
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results  in  the  following  expressions  for  the  Z- parameters  (Z  decreases 

with  depth) : 

ZC=Z  (7) 

Z  =  Z  +  K(S)  (8) 

Zb  =  Zc  -  K(S)  (9) 

Typical  values  of  K  are  1.47  and  1.30  for  temperature  increments  of  0.5 
and  1.0  degree  respectively. 

1 .   Non-Gaussian  Distributions 

The  technique  summarized  above  may  be  applied  to  non-Gaussian 
distributions  by  first  recalling  the  definition  of  Z   presented  at  the 
beginning  of  this  section.   As  noted,  the  original  data  must  be  re- 
structured to  define  a  distribution  function  Z(T)  -  the  ordered  array 
of  depths  read  at  a  constant  temperature  interval.   Once  this  is  done, 
Z   is  determined  by  straightforward  finite-difference  operations  on 
Z(T).   Z   corresponds  to  the  depth  at  which  the  second  finite  differ- 
ence passes  through  zero. 

With  Z  determined,  the  procedure  is  to  consider  separately 
the  segments  of  the  distribution  curve  lying  above  and  below  Z  .   For 

convenience,  let  the  functions  F  and  G  represent  the  two  segments  such 

< 
that  F  is  defined  for  Z  =  Z   and  conversely  for  G.   At  this  point  an 

artifice  is  introduced  whereby  each  of  the  functions  F  and  G  are  supplied 

with  a  "tail"  (to  be  graphic  if  not  exactly  precise)  fashioned  so  as  to 

transform  F  and  G  into  Gaussian  Thermoclines .   Standard  deviations 

(Sf  and  S  )  for  each  distribution  are  then  calculated  and  substituted 

into  equations  (8)  and  (9),  Sf  defining  Z  and  S   defining  Z  . 
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2.   Application  to  AXBT  Data 

In  the  analysis  of  the  AXBT  data,  the  temperature  distribution 
was  assumed  to  be  non-Gaussian  in  those  traces  with  one  major  thermo- 
cline.   In  the  case  of  well-defined  multiple  thermoclines ,  the  Z-param- 
eters  were  calculated  for  the  most  shallow,  assuming  a  Gaussian 
distribution.   Temperature  inversions  were  neglected. 

This  somewhat  cavalier  treatment  of  multiple  structure  does  not 
imply  such  features  were  considered  unimportant.   To  the  contrary, 
they  merit  a  much  more  rigorous  examination  than  the  scope  of  this 
thesis  permits. 

B.   LATITUDINAL  VARIATION  OF  THE  THERMOCLINE 

Figures  9-13  show  the  variation  of  the  top,  center  and  bottom  of  the 
thermocline  with  latitude  for  each  of  the  five  sections.   Again,  the 
seasonal  difference  in  gross  features  is  most  apparent.   Figures  9-12 
represent  conditions  at  the  end  of  the  heating  season  and  Fig.  13 
represents  conditions  after  about  two  months  of  cooling. 

The  dominant  feature  in  August  and  September  is  the  "dome"  in  all 
curves  in  the  mid-latitudes.   By  November,  this  dome  has  been  replaced 
by  a  distinct  "dip"  at  about  40  N  latitude.   The  explanation  for  these 
gross  features  in  the  absence  of  definitive  climatic  and  salinity  data 
is  speculative;  however,  the  features  are  not  inconsistent  with  the 
seasonal  mixing  processes  characteristic  of  the  subtropicand  subarctic 
regions  of  the  North  Pacific  (Tully,  1964).   In  the  subtropic  region, 
evaporation-driven  convective  mixing  is  predominant  throughout  the  year. 
In  the  subarctic  region,  where  precipitation  exceeds  evaporation,  the 
mixing  process  differs  in  the  heating  and  cooling  seasons.   During  the 


heating  season,  a  near-surface  halocline  created  by  excess  precipitation 
prevents  convective  mixing  so  that  the  mixed  layer  is  primarily  a 
function  of  wind  speed  (Tully,  1964).   As  the  water  cools,  however, 
temperature-driven  convection  occurs  despite  the  excess  precipitation. 
Keeping  these  mechanisms  in  mind,  it  appears  possible  that  the  dome 
present  in  Figs.  9-12  corresponds  to  a  transition  area  between  the  two 
regions  wherein  evaporation  and  precipitation  are  roughly  in  balance  and 
and  the  mixing  is  due  to  wind  alone.   By  November  (Fig.  13)  convective 
mixing  has  overcome  the  restraining  influence  of  the  seasonal  halocline 
north  of  40  N  latitude.   The  location  of  the  dip  in  Fig.  13  corresponds 
closely  with  the  southern  boundary  of  the  subarctic  region  defined  by 
Tully.   The  dip  suggests  that  the  intensity  of  the  gryral  circulations 
characteristic  of  the  subtropic  and  subarctic  regions  have  intensified 
to  form  a  convergence  zone  at  the  boundary,  thus  accelerating  the  down- 
ward mixing. 

Smaller-scale  features  are  more  difficult  to  discern  with  any 
certainty.   As  nearly  synoptic  as  the  sections  are,  the  fact  remains 
that  the  observations  for  a  given  transit  were  not  simultaneous;  thus, 
both  space  and  time  effects  are  present  in  the  data.   Further,  the 
spacing  between  observations  was  such  that  scales  on  the  order  of  a  few 
miles  are  nothing  more  than  noise  imposed  on  scales  on  the  order  of 
hundreds  of  miles  (i.e.  approximately  ten  times  the  station  spacing). 
Features  most  likely  to  evince  these  constraints  would  be  something  on 
the  order  of  large  thermal  fronts.   A  thermal  front  in  the  ocean  is 
loosely  defined  here  as  a  region  of  steep  horizontal  temperature  gradient; 
the  front  may  occur  at  the  surface  or  below.   A  more  quantitative 
definition  is  difficult  in  view  of  the  wide  range  of  intensities  assigned 
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to  thermal  fronts  in  the  literature.   Knauss  (1957),  for  example, 
describes  a  front  in  the  equatorial  Pacific  where  the  temperature 
changed  3  C  in  a  horizontal  distance  of  less  than  100m.   The  more 
general  case,  however,  is  usually  less  dramatic.   For  any  given  front, 
the  causitive  forces  are  not  always  apparent;  however,  the  presence  of 
a  front  is  usually  reflected  in  a  marked  change  in  the  spacial  distri- 
bution of  thermal  parameters. 

Excessive  smoothing  of  the  results  of  the  Z-parameter  calculations 
would  tend  to  mask  the  "discontinuities"  or  step-like  pattern  associ- 
ated with  oceanic  frontal  systems.   Therefore,  the  development  of  Figs. 
9-13  was  more  subjective  than  mathematical;  that  is,  envelopes  fitted 
about  points  rather  than  least-square  curves  were  examined.   Because  of 
this  dependence  on  personal  interpretation,  the  scatter  diagrams  from 
which  Figs.  9-13  were  developed  are  included  as  Appendix  A. 
1.   In  the  Heating  Season 

In  addition  to  the  dome  already  mentioned,  another  feature 
that  appears  regularly  during  the  heating  season  is  the  slight  upward 
trend  in  all  curves  just  north  of  50  N  latitude.   The  hazzards  of 
speculating  about  current  effects  based  solely  on  the  geostrophic 
component  are  appreciated;  however,  it  is  interesting  to  note  that  the 
upward  trend  corresponds  to  the  location  of  the  weak  westward  flow 
noted  in  the  velocity  patterns.   The  implication,  as  stated  earlier, 
is  that  the  westward  flow  results  from  a  convergence  of  the  Alaska 
stream  and  the  subarctic  current;  that  is,  the  Alaska  gyre  extends  at 
least  as  far  west  as  158  W  longitude. 

More  specific  details  of  thermocline  variation  and  the  relation 
to  surface  flow  during  the  heating  season  are  presented  below. 
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a.  19  August  (Fig.  9).   The  dome  structure  extends  between 
35-45  N.   The  vertical  extent  of  the  thermocline  is  smallest  in  this 
area  -  on  the  order  of  15m.   Fronts  appear  in  the  top  and  center  curves 
at  33  N  and  in  the  bottom  at  35  N.   The  top  and  center  frontal  zones 
correspond  to  the  shear  in  current  indicated  in  Fig.  3;  however,  other 
strong  shear  flows  have  no  apparent  effect  on  the  thermocline. 

b.  22  August  (Fig.  10).   The  upper  and  center  fronts  at 
33  N  remain,  although  the  center  has  deepened  slightly.   The  bottom 
front,  however,  has  deepened  about  18m  and  shifted  one  degree  to  the 
south.   A  current  shear  is  again  located  at  33  N.   An  additional  front 
has  developed  in  the  upper  curve  at  45  N;  in  the  same  time  period,  the 
surface  current  at  that  location  decreased  from  four  cm/sec  to  zero. 

c.  2  September  (Fig.  11).   The  upper  front  remains  at  45  N; 
all  other  frontal  zones  have  vanished.   The  dome  has  shifted  to  the 
south  about  five  degrees  and  appears  quite  smooth.   The  vertical  extent 
of  the  thermocline  is  still  smallest  in  the  dome;  however,  the  entire 
feature  has  deepened  about  four  meters.   There  is  no  apparent  correlation 
between  surface  currents  and  thermal  structure  anywhere  along  the  track. 

d.  4  September  (Fig.  12).   The  gross  features  along  the 
track  are  similar  to  preceeding  sections;  however,  the  change  in  detail 
developed  in  two  days  is  quite  pronounced.   The  dome  structure  has 
thickened  to  about  25m  and  a  "kink"  or  secondary  dome  has  developed  in 
the  south.   This  contrasts  to  the  smooth  appearance  of  the  dome  in 
previous  sections.   Fronts  have  formed  in  the  top  and  center  curves  at 
42  N;  the  center  front  is  quite  steep  and  corresponds  roughly  to  the 
sudden  increase  in  the  eastward  current  observed  in  Fig.  6.   These 
sudden  changes,  particularly  south  of  30  N,  lead  one  to  suspect 
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increased  atmospheric  effects  at  the  surface,  possibly  high  winds  or 
rain. 

2 .   In  the  Cooling  Season 

Figure  13  indicates  that  the  seasonal  thermocline  has  all  but 
disintegrated  and  is  sinking  to  the  level  of  the  permanent  thermocline. 
This  deterioration  is  accompanied  by  a  general  increase  in  eastward 
flow  between  25-50  N.   As  noted  earlier,  the  axis  of  the  dip  seems  to 
correspond  to  the  transition  area  between  the  subtropic  and  subarctic 
regions . 

A  large  smooth  "hump"  has  developed  in  all  curves  in  the  south. 
The  hump  seems  unrelated  to  the  dome  noted  during  August  and  September; 
however,  its  cause  is  not  clear.   It  may  be  the  result  of  a  Taylor 
column  effect  to  be  discussed  in  the  next  section. 

C.   IATITUDINAL  VARIATION  OF  ISOTHERMAL  LINES 

Figures  14-18  present  the  pattern  of  isotherms  along  the  five 
sections  together  with  the  traces  from  which  they  were  developed. 
Unfortunately,  the  scale  of  the  figures  eliminate  much  of  the  finer 
detail,  particularly  in  the  surface  layer.   However,  variations  noted 
in  the  thermocline  parameters  appear  to  be  consistent  with  the  near- 
surface  isotherm  patterns  when  examined  on  larger-scale  graphics. 

One  of  the  most  prominent  features  of  the  isotherm  contours  is 
the  abrupt,  frontal  character  of  the  8  -isotherm  -  rising  some  180m 
in  approximately  30  nautical  miles.   There  is  little  variation  in  the 
location  of  the  steepness  of  the  front  in  this  isotherm  below  about  60m. 
Neumann  and  Pierson  (1966),  among  others,  have  described  the  8  -isotherm 
as  the  boundary  between  the  cold  and  warm  spheres  in  the  ocean.   In  Figs. 
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14-18  it  seems  to  be  the  index  or  seasonal  activity  in  the  subarctic 
zone.   During  the  transition  from  the  heating  to  cooling  season  (Figs. 
17  and  18),  the  8  -isotherm  in  the  surface  layer  has  shifted  southward, 
intersecting  the  surface  at  about  48  N  in  November.   The  6  -isotherm, 
on  the  other  hand,  has  maintained  its  general  shape  but  deepened  about 
20m.   The  effect  of  the  steep  front  in  the  8  -isotherm  on  the  surface 
layer  does  not  reflect  in  the  figures;  however,  during  the  heating 
season,  the  front  seems  to  mark  the  northern  limit  of  the  dome  noted 
in  the  thermocline  patterns. 

The  10  -isotherm  generally  conforms  to  the  8  -isotherm.   In  over- 
all shape,  the  isotherms  south  of  (warmer  than)  10  C  are  more  horizontal 
from  10  C  and  colder,  the  isotherms  rise  abruptly  to  the  surface  layer. 
During  the  heating  season,  the  colder  isotherms  rise  to  between  20-40m 
and  then  run  roughly  parallel  to  the  surface.   The  pattern  of  the  8  - 
and  10  -isotherms  in  August  and  September  (Figs.  14-17),  particularly 
north  of  about  48  N,  suggests  the  lower  limit  of  an  estuarine  embay- 
ment  such  as  described  by  Tully  and  Barber  (1960). 

The  November  sections  (Figs.  13  and  18)  indicate  the  transition  area 
between  the  subarctic  and  subtropic  regions  to  be  between  41-45  N.   The 
wide  separation  between  the  10  and  12   isotherms  corresponds  to  the 
location  of  the  dip  in  Fig.  13.   Further,  the  shape  of  the  isotherms 
10  C  and  colder  is  noticeably  different  than  the  warmer  lines.   This 
suggests  the  10  -isotherm  as  the  southern  boundary  of  the  subarctic 
region. 

In  the  south,  the  limit  of  the  dome  noted  in  the  thermocline 
patterns  seems  to  correspond  to  the  18  -isotherm.   The  frontal  character 
of  this  isotherm  is  not  as  intense  as  the  8   line;  however,  the  southern 
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fronts  in  the  thermocline  curves  in  August  (Figs.  9  and  10)  coincide 
with  a  marked  rise  in  the  18  -isotherm  on  those  dates. 

The  most  interesting  feature  of  the  deeper  thermal  structure  is  the 
dome  or  "hump"  in  the  isotherms  between  23-26  N  latitude.   The  mechanism 
generating  this  feature  is  not  clear.   One  possibility  is  that  it  is 
simply  an  internal  wave  phenomenon  resulting  from  waves  either  reflected 
or  generated  by  the  Hawaiian  Islands.   Amphidromic  tide  waves  are  pro- 
gressive waves  that  can  be  reflected  by  obstacles  such  as  the  Hawaiian 
Islands.   There  is  general  agreement  that  the  semidiurnal  wave  approaches 
Hawaii  from  the  northeast.   However,  Wyrtki  and  Graefe  (1967)  point  out 
disagreement  as  to  the  direction  of  the  diurnal  wave  approach  to  Hawaii. 
The  point  is  that  the  progression  of  amphidromic  waves  across  the  North 
Pacific  is  not  certain,  implying  that  the  reflection  of  these  waves  is 
equally  uncertain.   Nevertheless,  the  possibility  remains  although  the 
consistency  of  the  dome  in  time  and  space  militates  strongly  against 
such  an  argument. 

A  second  possibility  is  that  the  feature  is  caused  by  a  Taylor- 
column  effect  such  as  described  by  Greenspan  (1968).   The  Taylor 
column  was  predicted  theoretically  by  J.  Proudman  in  1916  and  con- 
firmed in  the  laboratory  by  G.  I.  Taylor  in  1921.   The  essential 
feature  of  a  Taylor  column  is  a  fluid  system  in  a  state  of  solid  body 
rotation  in  which  a  fluid  motion  is  induced  relative  to  a  fixed  obstacle 
in  the  system.   To  illustrate,  consider  an  object  such  as  a  hocky  puck 
secured  to  the  bottom  of  a  filled  container  of  fluid  that  is  rotating 
on  a  turntable.   By  varying  the  rotation  of  the  turntable  slightly,  a 
slow  relative  motion  between  the  fluid  and  the  hockey  puck  is  estab- 
lished.  A  two-dimensional  flow  pattern  results  that  is  most  striking. 
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An  excellent    photograph   presented   by  Greenspan    (1968,    Fig.    1.2)    clearly 
shows   a    cylindrical   column   of    fluid   rigidly  attached   to  a    small    pro- 
tuberance   on   the    bottom   of  a   beaker.      Additional   photographs   demon- 
strate   that    flow  about    the    column    is   much    the    same   as    if   the    column 
were   an    impermeable   solid. 

Consider   now   the    rugged   bottom   topography   illustrated    in  Fig.    1   at 
approximately   the    same    latitudes    in  which   the    isotherm  humps   appear. 
The   urge    to    infer   a    relation  between   the    two    features    is    virtually 
irresistable   albeit    speculative. 

In  the  previous  section,  reference  was  made  to  a  possible  Taylor- 
column  effect  in  connection  with  the  appearance  in  southern  latitudes 
of  a  large  smooth  hump  in  the  thermocline  patterns  in  November  (Fig. 
13).  Recall  that  the  development  of  this  feature  was  coincident  with 
the  general  deterioration  of  the  thermocline.  The  implication  is  that 
if  a  Taylor-column  effect  does  exist,  it  is  bounded  at  the  top  by  the 
thermocline . 

Unfortunately,    there    is    little   evidence   documenting    the   existence 
of  Taylor   columns    in  non-homogeneous    fluid   systems    such  as    the    ocean; 
a    thorough    investigation   of   this    point    is    beyond    the    scope    of  this 
thesis.      However,    the    consistency   of   the    isotherm  humps    in    time    and 
space    is    sufficiently   compelling   to   suggest    that    the    possibility   of  a 
Taylor-column  effect   exists   and   merits    further   examination. 
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V.   THE  AXBT  AS  A  SCIENTIFIC  INSTRUMENT 

The  AXBT  was  developed  primarily  as  a  tactical  device  for  use  by 
antisubmarine  warfare  forces.   Its  purpose  is  to  provide  real-time 
intelligence  about  the  shape  of  the  temperature  versus  depth  curve  for 
use  as  input  to  various  range  prediction  techniques.   Hopefully,  the 
preceeding  sections  have  indicated  the  potential  of  the  AXBT  as  an  aid 
to  the  synoptic  collection  of  oceanographic  data.   It  remains  to 
determine  whether  or  not  the  data  so  collected  can  withstand  the  rigors 
of  scientific  analysis. 

Section  II  outlined  the  operation  of  the  AXBT.   In  this  section  the 
capability,  reliability  and  accuracy  of  the  instrument  will  be  dis- 
cussed.  Statistics  presented  in  this  section  were  determined  during 
an  operational  evaluation  of  the  AXBT  and  NADC  converter  by  the  U.  S. 
Navy  (COMOPTEVFOR,  1965). 

A.   CAPABILITY 

AXBT  were  launched  from  aircraft  at  altitudes  from  500  to  20,000 
ft.  and  groundspeeds  from  200-360  knots.   At  20,000  ft.,  however, 
complete  thermograms  were  received  from  only  33$  of  the  instruments 
dropped;  at  10,000  ft.  there  was  an  improvement  to  64$.   The  control- 
ing  factor  was  the  range  from  the  aircraft  to  the  instrument  at  the  end 
of  data  transmission.   Consequently,  it  was  determined  that  altitudes 
and  airspeeds  should  be  selected  such  that  the  aircraft  is  within  30 
miles  of  the  instrument  during  transmission  in  order  to  receive  a 
complete  thermogram. 
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The  NADC  converter  assembly  is  sufficiently  small  and  light  for  use 
in  helicopters.   An  added  advantage  of  the  system  is  that  rotary-wing 
aircraft  need  not  hover  during  data  collection. 

B.  RELIABILITY 

Seventy-six  percent  of  137  AXBT  expended  during  the  Navy  tests  were 
considered  to  have  provided  satisfactory  thermograms.   Most  mal- 
functions were  connected  with  the  audio  oscillator  controlled  by  the 
temperature  probe;  common  failures  were  abrupt  frequency  shifts,  early 
termination  of  the  signal,  and  weak  signal. 

Most  malfunctions  are  easily  recognizable  by  trained  operators.   One 
exception  is  a  hung  probe  in  isothermal  water;  both  conditions  produce 
similar  traces. 

C .  ACCURACY 

Accuracy  of  the  AXBT,  compared  to  the  mechanical  BT,  was  determined 
during  five  separate  tests  with  a  total  of  84  AXBT.   The  airborne 
instruments  were  dropped  within  500yd  of  an  oceanographic  research 
ship.   Mechanical  thermograms  were  obtained  by  the  surface  ship 
simultaneously  with  the  airborne  where  practicable,  but  in  no  case 
was  the  comparison  BT  more  than  one-half  hour  old.   The  airborne 
thermograms  were  those  displayed  by  the  NADC  converter  assembly. 

A  total  of  1,141  data  points  from  the  airborne  thermograms  were 
compared  to  the  reference  data.   The  mean  algebraic  error  was  0.032  F 
with  a  standard  deviation  of  1.039  F.   The  mean  absolute  error  was 
0.581  F  with  a  standard  deviation  of  0.861  F.   These  figures  reflect 
contributions  from  the  AXBT  and  the  NADC  converter  as  well  as  possible 
errors  induced  by  time  and  space  separation  between  the  AXBT  and  the 
reference.   The  magnitude  of  each  contribution  is  not  known. 
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In  another  test,  a  comparison  of  32  distinctive  breakpoints  noted 
on  airborne  thermograms  with  the  same  breakpoints  from  reference 
thermograms  showed  that  88$  agreed  within  20  ft.  in  depth. 

D.   APPLICATION 

As  a  result  of  the  tests  described  above,  the  AXBT  was  recommended 
for  acceptance  by  the  Navy  in  late  1965.   The  accuracy  of  the  instrument 
as  compared  to  the  mechanical  BT  was  considered  satisfactory  for  the 
tactical  nature  of  its  intended  use.   Implicit  in  this  acceptance  by 
the  Navy  is  that  the  AXBT  enables  the  determination  of  acoustic 
parameters  such  as  layer  depth  with  sufficient  accuracy  and  speed  to 
maintain  hot  pursuit  of  submerged  contacts. 

The  value  of  the  instrument  as  a  scientific  tool  depends  greatly 
on  the  nature  of  the  problem  under  investigation.   The  instrument  appears 
to  provide  reasonably  consistent  and  reliable  information  about  relative 
structure  (for  which  it  was  designed);  however,  the  mean  absolute  error 
reported  above  (0.581  F  or  about  0.32  C)  poses  a  serious  limitation 
to  studies  requiring  precise  values  of  temperature.   For  example,  the 
direction  of  geostrophic  velocities  calculated  by  both  the  dynamic  and 
T-S  gradient  methods  showed  considerable  variability  when  calculated 
for  stations  spaced  about  50km  apart;  the  variability  decreased  with 
increased  separation.   The  implication  is  that  the  increased  temperature 
differential  lessened  the  effect  of  the  error.   This  points  to  the  real 
value  of  the  AXBT  -  large-scale  coverage.   The  AXBT  and  its  shipborne 
counterpart,  the  XBT,  represent  the  first  true  improvements  in  thermo- 
graph instrumentation  since  the  introduction  of  the  mechanical  BT  in 
1937.   Prior  to  the  development  of  the  expendable  BT,  instruments 


57 


employing  electronic  circuits  rarely  were  able  to  utilize  the  best  in 
electronic  circuitry  either  due  to  prohibitive  cost,  bulk,  or  size  of 
components.   Dramatic  advancements  in  the  field  of  microelectronics 
generated  by  the  space  program  served  to  ameliorate  many  of  these 
difficulties . 

Microelectronics  appears  to  offer  the  best  means  for  improving  the 
AXBT.   Assuming  a  sufficiently  wide  acceptance  of  the  instrument  to 
provide  a  broad  base  for  amortizing  the  cost,  the  AXBT  could  be  greatly 
improved  by  incorporating  an  additional  thermal  sensor  at  the  surface. 
This  should  provide  a  more  accurate  sea  surface  reference  temperature 
than  is  presently  recorded  by  the  temperature  probe  during  the  75 
seconds  it  remains  at  the  surface. 

The  most  immediate  obstacles  to  be  overcome  before  wide  acceptance 
of  expendable  instruments  by  the  scientific  community  are  more  subtle; 
they  involve  people  and  psychological  reactions  to  the  philosophy  and 
concepts  represented  by  the  AXBT.   When  using  the  mechanical  BT ,  tnean- 
ingful  results  are  greatly  dependent  on  careful  handling  and  careful 
calibration;  calibration  being  affected  as  often  as  practicable.   The 
production  techniques  used  in  the  manufacture  of  the  AXBT,  however, 
are  such  that  the  instrument  must  be  calibrated  in  the  factory  with 
no  provision  for  field  calibration.   One  can  appreciate,  then,  the 
reluctance  of  the  oceanographer  to  trust  an  instrument  which  he  does 
not  calibrate  himself. 

Confidence  in  instruments  such  as  the  AXBT  and  XBT  can  only  come 
with  time.   Improved  electronic  production  techniques  involving  care- 
ful application  of  quality  control  in  reliability  engineering  should 
make  it  possible  to  assign  a  confidence  figure  to  the  performance  of 
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each  instrument  (Snodgrass ,  1966).   The  goal,  of  course,  is  to  develop 
confidence  in  the  instrument  to  the  point  where  credibility  can  be 
assigned  to  temperature  idiosyncracies  observed  even  in  a  single  trace 
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VI.      CONCLUSION 

A.  SUMMARY 

The  motive  behind  the  approach  to  oceanic  variability  pursued  in 
this  thesis  has  been  to  determine  the  amount  and  reliability  of  infor- 
mation that  could  be  determined  solely  from  temperature  data.   Latitudi- 
nal distribution  of  thermocline  depths  and  isothermal  lines  were  selected 
for  examination  in  order  to  demonstrate  thermal  structure  variability  in 
the  most  graphic  form  possible.   Various  techniques  were  attempted  to 
determine  thermocline  depths;  the  quasi-statistical  approach  of  the 
Gaussian  thermocline  produced  the  most  satisfactory  scatter  diagrams 
and  was 3  therefore,  selected  for  presentation  here.   A  cursory  examin- 
ation of  geostrophic  velocity  fields  was  possible  because  of  the  recent 
introduction  of  the  T-S  gradient  method. 

B.  CONCLUSIONS 

The  consistency  of  the  data  evidenced  in  the  five  space  sections 
developed  from  airborne  bathythermograms  appears  quite  good.   Based 
on  the  examination  of  these  sections  presented  earlier,  the  following 
conclusions  are  offered: 

1.  The  T-S  gradient  method  is  a  viable  technique  for  deter- 
mining geostrophic  velocity.  It  is  computationally  less  complicated 
than  the  dynamic  method  and  provides  comparable  results. 

2.  The  definition  and  calculation  of  thermocline  parameters 
using  the  notion  of  a  Gaussian  thermocline  is  a  more  rigorous  and 
scientific  approach  to  the  problem  than  the  plethora  of  arbitrary 
gradient  methods  that  abound  in  the  literature.   The  technique  is  not 
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restricted  to  thermal  parameters  or  to  Gaussian  distributions.   The 
variation  with  depth  of  any  oceanographic  parameter  (such  as  salinity, 
dissolved  oxygen,  etc.)  is  adaptable  to  a  similar  analysis  once  a 
reasonable  statistical  representation  is  selected,  for  example  the 
chi-square  or  student's  t  distributions. 

3.  The  variation  in  depth  to  the  thermocline  during  the  heating 
season  is  characterized  by  a  distinct  dome-shaped  distribution  with 
latitude . 

4.  A  disturbance  in  the  water  column  exists  between  23  N  and 
26  N  latitude  as  evidenced  by  the  persistent  "hump"  in  the  distribution 
of  isothermal  lines  with  latitude.   The  presence  of  a  Taylor-column 
effect  is  possible. 

5.  Thermal  fronts  tend  to  form  between  33-35  N  and  43-45  N 
latitude  during  the  heating  season. 

6.  The  AXBT  has  great  potential  as  a  scientific  instrument  and 
represents  a  significant  step  toward  synoptic  oceanographic  data 
collection. 

C.   RECOMMENDATIONS 

The  PARKA  data,  considered  in  total,  offer  sufficient  flexibility 
in  time,  space  and  scale  to  accommodate  a  variety  of  studies  too 
numerous  to  itemize  here.   Further,  the  treatment  of  the  AXBT  data  has 
been  far  from  exhaustive.   However,  it  is  left  to  the  reader  to  specu- 
late on  the  full  import  of  so  large  a  body  of  data.   Recommendations 
for  further  study  enumerated  below  are  restricted  to  those  points 
specifically  related  to  the  development  of  the  data  as  presented  in 
this  thesis.   Specifically: 
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1.  The  possible  existence  of  a  Taylor-column  effect  in  the 
waters  just  north  of  Hawaii  suggests  the  need  for  a  theoretical  study 
of  a  simple  two-layer  model  to  investigate  the  behavior  of  a  Taylor 
column  in  a  stratified  fluid.   As  a  minimum,  the  study  should  address 
questions  such  as  : 

Does  the  column  have  any  effect  on  the  upper  layer? 

Can  the  column  withstand  an  induced  current  profile? 

Are  there  critical  dimensions  for  the  obstacle? 
Eventually,  the  model  should  be  expanded  to  include  light  stratification 
of  the  lower  layer. 

2.  The  results  of  geostrophic  velocity  calculations  support  the 
contention  by  Reed  and  Laird  (1966)  that  station  spacing  is  a  critical 
parameter.   Further,  it  is  reasonable  to  expect  optimum  spacing  to  be 

a  function  of  measurement  error  in  determining  the  independent  vari- 
ables in  the  velocity  equation.   Accordingly,  it  is  recommended  that 
a  study  similar  to  that  presented  by  Reed  and  Laird  be  conducted  in 
order  to  establish  optimum  spacing  criteria  for  use  with  the  T-S 
gradient  method. 

3.  It  is  recommended  that  statistical  curve-fitting  techniques 
similar  to  Boston's  Gaussian  thermocline  be  developed  for  the  analysis 
of  other  oceanographic  variables.   In  order  to  increase  computor  storage 
capacity,  particular  emphasis  should  be  placed  on  determining  the  mini- 
mum number  of  statistical  parameters  required  to  reproduce  a  given  dis- 
tribution.  Implicit  in  this  recommendation  is  a  re-examination  of 
digitizing  procedures  wherein  values  of  a  selected  parameter  are 
recorded  at  predetermined  depths  ;  the  reverse  process  described  in 
Section  IV  may  be  more  appropriate. 
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4.   Finally,  it  is  recommended  that  a  design  study  be  conducted 
to  determine  the  feasibility  of  providing  the  AXBT  with  an  additional 
thermal  sensor  at  the  surface  in  order  to  increase  the  reliability  of 
the  sea-surface  reference  temperature. 
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APPENDIX  A 

THE RMOC LINE  SCATTER  DIAGRAMS 

Figures  19-23,  on  the  following  pages,  present  the  scatter  diagrams 
from  which  figures  9  to  13  were  developed.   The  scatter  diagrams  are 
presented  solely  to  afford  the  reader  an  opportunity  to  evaluate  the 
author's  judgement  in  the  subjective  development  of  the  curves 
describing  the  latitudinal  variation  of  the  top,  center  and  bottom 
of  the  thermocline. 

The  depth  scale  in  Figs.  19-23  is  arbitrary.   The  scatter  of  points 
for  the  three  z-parameters  have  been  separated  in  depth  for  clarity; 
refer  to  figures  9-13  for  approximate  depths. 
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